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ABSTRACT 



A study of a multimode third order instrument servomechanism is 
presented. Particular emphasis is placed on the effects of the initial 
condition of a compensating network, switched in as the system enters 
a linear operating mode, upon the behaviour of the output. The optimum 
values for the initial condition as a function of velocity and accelera- 
tion of the servomotor is derived, using the integral squared error as 
a criteria. The results are confirmed with a computer simulation. 
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List of Symbols and Abbreviations 



x(t) 

x(t) 



= Damping factor of a complex pole second order system 
= Natural frequency of the system 
= Time constant of a system, viz.: = RC 

= Bode gain of type one system, or Velocity constant 

= dx(t) 
dt 

d2x(t) 

, . 2 
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I. INTRODUCTION 



Multimode control systems are those which are operated under two 
or more controlling laws which are in turn switched on or off according 
to some pre-established rule. Here, a particular case is studied, in 
which an instrument-servo is used to position a device in two main 
steps: 1,- First the motor is driven open loop at full forward speed, 

and then slowed down by reversing its voltage — both without feedback; 
2,- Finally, within a given error limit, a feedback loop is closed, 
zeroing the error according to the linear-system response. 

It is the second of these steps that is studied with greater detail 
as it is here where the optimization technique of this work applies. 

In general, in this type of system a third order type one instrument 
servo is used which is unstable under the desired gain conditions. 
Therefore, a compensating network has to be added to the forward path 
of the loop which is a lead network in this case. In practice, the 
initial condition on the capacitor of this lead network is usually made 
equal to zero, as it is easier to have a short circuit across the capac- 
itor than to provide an extra voltage source for this initial condition. 
No investigation has been found in the literature as to the effect of 
this initial condition on the time response of this type of control, and 
it is the purpose of this study to determine its effect on the servo 
performance. 

To accomplish this, a typical third order instrument servo of the 
linear actuator type has been modeled and the whole problem simulated 
with the aid of the System/ 360 Continuous System Modeling Program of 
IBM ( S/360 CSMP ) in the IBM 360/67 Computer. 
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The results obtained from this special case simulation will be 
discussed in the light of producing a more general criteria for similar 
cases . 

The performance criteria to be used will be that of the Integral 
Squared Error, since this function provides the necessary information 
on minimum error which is desired here. Minimum settling time will be 
observed from the data obtained from the computer simulation. 
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II. BASIC CONCEPTS ON DISCONTINUOUS SYSTEMS 



Consider first a general block diagram of a discontinuous control 
system as shown on Fig. 2.1. The general control is achieved by the 
control computer which can be an analog or digital computer, an ampli- 
fier with relay, or, a simple relay circuit. The controlling signals 
may come from any point in the system, the input and output being the 
most frequently used. The control computer uses its input signals to 
provide the desired outputs which control the parameters of the plant, 
the controller and amplifier, or the feedback signals. 

From the analysis of the input signals, the input signal analyzer 
provides the required signals for the control computer. This analyzer 
may consist of a differentiator or integrator network, a sampler, a 
potentiometer or a function generator. 

The controller and amplifier network may consist of an amplifier 
with relay with R-C network, or a low power amplifier circuit which 
provides a compensated signal for the plant drive. The control computer 
can control the gain of the amplifier and the compensating circuit. 

The plant is the object to be controlled by the rest of the devices. 
In electrical control systems it may consist of a motor driving a load, 
its shaft being designed to follow the input signal. The plant must also 
be designed to meet the specifications of the system, which in the given 
case would be settling time, overshoot, transient response, etc. The 
control computer can vary the gain, the damping and/or the inertia of 
this plant. 

The feedback block consists of one or many signals which are a 
function of the outputs. The design will consider the number of signals 
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Fig. 2 ..]. General Block Diagram For a Discontinuous System 



to be used, and the control computer can control the polarity, instant 
and period for them to be applied, as well as their gains and character. 

From this general structure, many special cases can be derived by 
placing limitations or specifications on the blocks indicated above. 

By using the controller and amplifier block and the plant, various 
control methods have been developed. The basic theory along these lines 
is to design the controller to make the switch operate along a curve in 
the phase plane, or a surface in the phase space that is called the 
optimum switching curve or hypersurface (Ref. [1]). 

For the continuous operating mode of the system, a compensating 
network is used, and if blocks 2 and 3 are used, the compensation is the 
cascaded type. This form of compensation seems predominant in this 
field, although other types of compensation could be used. The theory 
of cascaded lead - lag networks has been largely discussed and developed, 
the Bode diagram being used with great advantages in the design. 

The reason for using discontinuous systems arises from the need of 
applying a limited amount of power to accomplish traveling in minimum 
time. These types of systems have been proven to be the best form of 
control even under economic considerations. 

Discontinuous systems are based on the theory that it is best to 
use maximum drive at the beginning of the trajectory, and then brake 
the system at a proper point. Thus, full voltage is applied at the 
beginning, and full reverse at the switching point. With the disconti- 
nuity, the character of the system is entirely different when the state 
point is at different sides of the hypersurface. The decaying part of 
the trajectory should be made, as close to the E vs E plane as possible 
Ref. [2]. 
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Comparing the continuous system with the discontinuous, the phase 
« 

plane plot of E vs E serves best to give a better insight of the differ- 
ence. Fig. 2.2 shows a typical phase plane of a discontinuous system, 
whilst Fig. 2.3 presents that of a second order continuous system. 

The continuous system is normally designed to have some overshoot 
in the trajectory, because an operation in the overdamped or even 
critical condition provides a slow response. In the discontinuous 
system the acceleration trajectory can be made unstable, and the die 
out trajactory a straight line near the E axis in the phase plane 
(fast eigenvector). 

This general discussion has only concerned itself with possible ways 
to use discontinuous theory to improve the response character of systems. 
The general switching theory and response character will not be discussed 
here in more detail. The main concern in later chapters will be the 
effect of entering a final phase which will be continuous; the use of 
a lead type conpensating network, and the influence of the initial 
conditions, in particular that of the compensating network. 
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Ftg. 2,2 Typical E vs E Plot for a Discontinuous Operating System 
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Ftg .. 2.3 Typical E vs E Plot for a Second Order System 
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HI. PLANNING THE EXPERIMENT 



A. ANALYTIC DISCUSSION 

In order to study a case of how the initial conditions on the lead 
network would affect the performance of a system, the problem was simu- 
lated using the characteristics of a linear actuator, such as the ones 
used in disc memories for computers to move the arm which carries the 
magnetic head. 

The system is supposed to move different distances, the only known 
variables being the initial error and the actual position of the head 
at all times. Minimum time and final error are of the essence, thus 
leading to a discontinuous bang-bang type of operation. The final part 
of the trajectory, i.e.: when the head reaches the edge of the track 
to which it was ordered, is accomplished in a closed loop mode. An 
extra gain is introduced to the system, and thus a filter has to be 
added in cascade in order to achieve stability, the filter being a lead 

network. 

Once a switching is determined, in order to have a real starting 
condition of velocity and acceleration of the motor as it enters the 
closed loop mode, a computer simulation using the S-CSMP/360 was 
implemented to study this final phase. 

1 . Equations of the Physical System 

The known characteristics of the motor in practical MKS units 

were: 

. . _ ,r r Newton 

Force constant = K m - Ib.b 

Motor inductance = L = .01 [Henry] 

a 
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Motor resistance = R = 5 [Ohms] 

a 

Mass of motor = M = 0.45 [Kg.] 

Newton -sec 

Viscous damping = f = 17.5 Metir 

Volts -sec 

Back EMF constant = K b = 15.6 — Meter 

Armature voltage = v fl [Volts] 

Mx + fx = K m i a 



(3.1) 

(3.2) 



which, using Laplace transforms and substituting for i a , gives 



V Cs) = 



Cs 2 M + sf) 



R s L o 
a + a (s 2 m + sf) + sK 

Sn K m 



X(s ) 



or 



G(s) 



= H * 1 = 

V(s) 



1 



[ s 



L M 
2 a 



m 



MR. L f fR 

mm m 



(3.3) 



which is the transfer function for the motor with voltage source feed. 

Similarly, the study could be made by considering a current 
source feed, in which case the transfer function would be 

X(s) = K m (3.4) 

I(s) s(sM + f) 

Going back to the voltage feed transfer function, and substi- 
tuting the given values of the motor characteristics, the open loop 
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gai n is : 



G(s) = 



3432.3 



s[s 2 + 538.6s + 72741] 



(3.5) 



From this expression, and by simple inspection, it is known 



that 



269.7 



' rad 
sec 



5 = 0.9985 



for the quadratic factor in the denominator. 

Thus the servomechanism is a type one, third order system, with 
two complex poles almost critically damped. The Bode diagram for such 
a motor is shown in Fig. 3.1. 

From the Bode diagram, it is seen that the K y of the motor is of 
the order of 0.047, which will represent an exceedingly slow system for 
the final closed loop phase of the operation, as compared with the desired 

settling time. 

2. Filter Design 

By further inspection of the Bode diagram, it is determined that 
an amplifier which should provide an extra gain of 55,000 ought to be 
used. This new characteristic function gives an unstable closed loop 
system, which in turn leads to the necessity of adding the lead network, 

A general and most simple form of lead network is represented by 



G c (s) - 



Ct(sT + 1) 
asx + 1 



(3.6) 



and its physical implementation is shown in Fig. 3.2, where a is the 
relationship between the corner frequencies of the filter as referred 
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Ft g - 3.2 Lead Filter Diagram 




Ft g ^ 3.3 Block Diagram of Implemented System 
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to the Bode plot, i .e. : 



a = 



R 



1 



R + R 



1 



In order to obtain the maximum phase correction with the filter, 
an a of 0,1 is used [3], and by trial and error process the system 
is compensated to provide the Bode plot shown in Fig, 3.4, the filter 
having the following characteristics 
a = .1 

t = 0,004 sec 



Thus 



G c (s) = 



0.1(.004s + 1) 
.004s + 1 



(3.7) 



The compensated system has now a transfer function 



G(s) 



5500(0. 004s + 1) 3432.3 

0.0004s +1 s(s 2 + 538.6s + 72741) 



(3.8) 



It is a stable system with a phase margin of 52 degrees, its 
root locus plot has been presented in Fig. 3.5. 

B. OBSERVATION OF THE EFFECT OF THE INITIAL CONDITIONS 

With the help of the flexibility of CSMP, several values of initial 
conditions for the filter can be implemented in one computer run, thus 
allowing a gross estimation of the direction of the variation of the 
performance with the said initial condition. 

The performance criterion used is that of the integral squared error 
in the linear mode. The optimum value of v^ is calculated to minimize 
this index. Since there are no weighting factors to consider in this 
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Fig, 3.5 Root-Locus Plot for the Compensated System 

xscale = 500 units/tnch 
yscale = 5 units/inch 



case, the problem is a very simple one, and the general expression for 
the optimal control 




(3.9) 



where t Q and t^ are the initial and final time; h and g are scalar 
functions; x(t) is the trajectory and u_(t) the control. £Ref. 4] 
reduces to simply 



and is given by the transfer function and solved by the program itself. 

Although flexible, CSMP places a number of restrictions on the use 
of subprograms and subroutines, therefore making it quite difficult to 
use the normal and known procedures for minimization and optimization. 
This leads to a rather crude method of determining the best suited 
initial condition for the filter, namely that of analyzing an ample 
range of values which show how the response of the system, and in 
particular the integral squared error, vary with the variable in 
question. Once a relative and feasible range of values for the initial 
condition is found, the best fit can be determined without much effort. 

1 , Algebraic Development 

The solution for x(t) will now be pursued using the state 
variable form approach. 



tf 

0 = / E-E dt 

0 



(3.10) 



in which 



E = g(x(t),u(t),t) 



(3.11) 



From the motor equations (3,1) and 3.2), it follows that: 



x(t) = -ax(t) - bx(t) + cv (t) 

a 



(3.12) 
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K K h fR. 

where a = pjj- — + = 72741 



b = 



I + r- = 538.6 
M L a 



c - mt = 3432 * 3 
6 



Also, from the filter circuit in Fig. 3.2, 
i(t) = i R (t) + i c (t) 



e ft) v c (t) 



+ C v c (t) 



e 0 (t) = e i Ct) + v c (t) 



Thus 



(3.13) 



(3.14) 

(3.15) 



vl_ J_ e.(t) 

v c (t) = . ( RC R,C ) v c (t) + — 

and since RC = T 
and R = 9R-| 

v c (t) = V c (t) 

This result holds under the assumption 



(3.16) 



(3.17) 

that the motor does not 



load the filter circuit, which can be achieved in practice with the 
use of an isolating amplifier. With this in mind, the analysis can be 
furthermore pursued to solve" -for the relationship between the output 
x(t) and the initial condition on the capacitor v^(0). 



Since the system is closed loop with zero input, the general 
expression 

E = R - C 



reduces to 



E = - C 



and therefore 

e.(t) = Ke(t) = - K. x(t) 

where K is the gain of the forward feed amplifier. 
From (3.17) and (3.18) 



(3.18) 



v c (t) = . ia v c ct) - ^ x(t> 



(3.19) 



Similarly, taking equations (3,12) and 3.15), and considering 

that 

e 0 (t) = v a (t) 

then: 

x(t) = -cKx(t) * ax(t) - bx(t) ~ cv^t) (3,20) 



Expressing equations 3.19 and 3.20 in their state 
variable form, and letting 
x 1 - x(t) 
x 2 = x(t) 
x 3 = x(t) 
x 4 = v c (t) 

then 

x(t) = Ax(t) + Bu(t) (3,21) 
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where 



A = 



0 

0 

*cK 

-9K 

T 



1 

0 

-a 



0 

1 

-b 

0 



0 

0 



- 10 . 

T 



and u(t) = 0 



Using Laplace transforms, the solution to equation (3.21) is 



-1 



x(t) = X [si - A]~ 1 x(0) 



from which, the desired solution is the expression for x-j = x(t) 
that is 



x i tt) =x 






1 



10 w _2. 



s(s 2 +bs+a)(s + — ) + cK(s + — ) 

T T 



— (s + — )(s +bs+a)x, 
+ I) x I 



(0) 



+ (s+b)(s + ~)x 2 (0) + (s + *^)x 3 (0) + cx 4 (0) 



xCt) -i 



-1 



_ (s+m)(s+n)[(s+p) 2 + q 2 ] 



(s + — ) (s 2 +bs+a)x(0) 

T 



+ (s+b)(s + — )x(0) + Cs + — )xC0) + cv r (0) 

T T 



where 



m 

n 

P 

q 



248,019 
2533.1968 
128,696 - 
241 ,98 
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Thus 



x(t) = Ea-jX(O) + a 2 x(0) + a 3 x(0) + a 4 v c (0)]e 



-mt 



+ [b ] x(0) + b 2 x(0) + b 3 x(0) + b 4 v c (0)]e 



-nt 



+ Ic^xCO) + c 2 x(0) + c 3 x(0) + c 4 v c (0)]e"^ lt (cos qt - p sin qt) 



+ [d^x(O) + d 2 x(0) + d 3 x(0) + d 4 Vj,(0)]e~P t sin qt (3.22) 



wi th 



a-j = 9.06x10 
a 2 = 3.9x10' 
a- = 1 .35x10 



-3 



-6 



■2.06x10 



-4 



b ] = 1.27x10 



-2 



bo = 



-4.96x10 v 
-9 



bo = 2.49x10 



b„ = 2.57x10 



-7 



= 9.78x10 



-1 



Co = -3.93x10 



-3 



Co " 



3.54x10 
-4 



-6 



Ca = 2.06x10 



d-j = 1.18 

do = 3.93x10 



-3 



d 3 = -5.01x10 
d. = 2.36xl0~ 4 



Optimization of v c (0) 

The general expression for the Integral Squared Error is 



ISE = / f (R-C) 2 dt 



(3.23) 



that in this particular case, with R = 0 and C = x , reduces to 



t f 2 - 

ISE = / x (t) dt 

0 
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To obtain a relative minimum of this function with respect to 
v^-(O), its partial derivative with respect to the variable is made equal 
to zero, therefore 



alSE _ l f 

t^oT ~ 0 / 




0 



(3.24) 



Thus, using equations (3.22) and (3.23) 



9ISE 



-mt 



-nt 



' avpto T = 0 / x t t )( a 4 e ""' 1 ' + V"*' + c 4 e " pt ( cos qt - p q sin qt) 
+ d^e^^sin qt) dt = 0 



If t^ is taken large enough, say of the order of the settling 
time, the equation for v^(0) can be greatly simplified, since the expo- 
nential terms tend to zero. Hence the integration gives: 

v c (0) = 25865. 3(a-jX(0) + a 2 x(0) + z 3 ^(0)) 

-14363.1 (b,x(0) + b 2 x(0) + b 3 ^(0) ) 

-29476. 3(c.|X(0) + c 2 x(0) + c 3 #(0)) 

Taking typical values of initial conditions as those obtained 

-5 

from the computer simulation, viz.: x(0) = -6.35x10 

x(0) = .27 
k0) = -190 

the optimum value of v^(0) is 79.3 volts, which corresponds to an 
initial condition of -79.3 volts. It will be shown later that this 
result agrees with the experimental values obtained from the computer. 
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To show that the inflection point of ISE vs v c (0) is in fact 

a minimum, the second derivative of the function is determined 

■ a - 2 . I ^ E - ^ = 2 (a.e"" 11 + b 4 e~ nt + c 4 e' pt (cos qt - ^sin qt) 

av c 2 (0) 44. h 

+ d 4 e" pt sin qt) dt 

= 4x1 Cf 2 

A similar analysis could have been made to determine the optimum 
value of any other of the initial conditions or a combination of them. 

The next section will be dedicated to analyze how the whole 
system is influenced by using different values of capacitor voltage. 
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IV. COMPUTER RESULTS 



A. DETERMINATION OF INITIAL CONDITIONS 

In order to have realistic valued initial conditions of the system 
when it enters the linear operational mode, a computer simulation of 
the bang-bang operation phase was accomplished. A convenient switching 
law was determined, such that for large enough errors, viz.: ten or 

more tracks, the motor was taken to around zero error with the appropri- 
ated final values of velocity and acceleration that would not produce a 
positive or negative overshoot when the loop was closed. 

It is important to note here that if the arm override enters a 
neighboring track, a new control law is encountered, and thus the system 
is driven completely out and away from its destination track. Therefore, 
it is imperative that once the head has entered the desired track zone 
it does not leave it, and if it does, it is only by a very small and 
insignificant distance (one mil or so), so that the inertia of the 
system would take care of the override, not allowing the head to be 
driven off. 

With this in mind, the results obtained from the discontinuous 

operation mode were tried as initial conditions on the closed loop mode. 

_2 

The width of the track considered was of 5 mils or 12.7x10 mm. 

Once a set of velocity and acceleration initial conditions were 
proved satisfactory to keep the head within the destination track, a 
rounded variation about these values was tried out in order to determine 
the permissible amount of error in these parameters in entering the 
linear zone. It should be remembered that there is no accurate means 
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to determine velocity or acceleration, hence the need of having a 
possible wide range of these values at the entering point, and a 
control law that is based on position alone. 

A computer program for the discontinuous phase is presented in 
the appendix, the typical phase plane plots for this mode being shown 
on Fig. 4.1 . 

B. TYPICAL PHASE PLANE RESULTS OF CONTINUOUS SYSTEM 

The computer program used in this part of the experiment is further 
explained in the appendix. More details on the programming method can 
be found in Ref. Jsj 

The simulation proved satisfactory to show the way the response of 
the system varies with the initial condition on the filter capacitor. 

By varying this initial condition over a wide range of values, the 
integral squared error as a function of the capacitor voltage v/as 
plotted, and the resulting curves are shown on Fig. 4.2. It can be 
seen from this figure that the integral squared error does minimize 
within the range of practicable voltages for the capacitor. Also, from 
the phase plane plots shown on Fig. 4.3, and the transient response on 
Fig. 4.4, the effect of the initial condition proves to be of great 
importance in obtaining the solution to the problem at hand. A zero 
initial condition turned out to be far from ideal in this case, and the 
voltage derived in the analytic discussion showed to be the best 
solution in practice. 

To have a better idea of the transient response of the system, a 
computer run for a unit step input and zero initial conditions was made 
the position vs. time curve being shown on Fig. 4.5. By a comparison 



30 




Fig. 4,1 Phase plane plots for discontinuous operating mode of system. 

The three curves show different values of positioning orders; 
the largest used was 400 tracks, the smallest 20 tracks. 
Switching law used is a linear one, namely E = INITIAL ERROR/7, 
x scale = .01 units/inch 
y scale = 0.2 units/inch 
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Fig. 4.2 Plot of integral squared error vs. filter initial condition 
for different pairs of velocity and o acceleration values. 
Note^that the optimum value for v c (°0) with the given x(t) 
and x(t) is around -60 volts. u 
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Ftg, 4,3 E vs E for x(0)=-160 m/sec 2 and x(0)=.12 m/sec. 
Curve (a) shows response for v^.(0)=-90 volts, 

(b) for v c (0)=-85 volts, and 

(c) for v q( 0)=-80 volts. 

x scale = ,02 m per inch 
y scale = 5x10”^ m/seg per inch 
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and acceleration, vq( 0)=0 in this case. 

x scale = .01 sec per inch 
y scale = 2x.0“ m per inch 
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Fig, 4,5 Step Input Response of the System. 

x scale = .01 sec per inch 
y scale = .2 m per inch 
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of this curve with the ones obtained as solutions to the discontinuous 
problem, it can be seen that the initial conditions on velocity and 
acceleration greatly affect the behaviour of the system, making its 
stability quite critical depending on their values. To control this 
situation, the initial condition on the filter capacitor has to be set 
to an appropriate value which is dependent on the velocity and accelera- 
tion of the motor at the point of entry to the linear mode. Therefore, 
the best possible solution would be one that used a different value of 
capacitor voltage for each set of initial conditions at this point. 

Again, this cannot be done because of the lack of information on these 
values, thus leading to the use of a best fit voltage that will apply 
to a whole set of parameters. 

It has been the intention to show here, with the help of the phase 
plane plots, of how the error varies with the initial condition on the 
capacitor. The various computer runs showed that, although this voltage 
could prevent a negative overshoot of the error at the beginning of the 
trajectory, it produced a positive one later on if not chosen properly, 
thus providing an undesired solution. In this respect, the integral 
squared error was very efficient in detecting these effects, and 
although with the specified initial conditions of acceleration and 
velocity in some cases the system did overshoot and left the bounded 
region, using the initial condition on the filter as indicated by the 
minimized integral squared error., the reader head was maintained within 
the specified limits. Fig. 4.6 shows the difference of the trajectories 
and the overshoot caused by the use of a wrong voltage, and how it is 
remedied with an adequate one-. The values of the error when overriding 
were of the order of 0.7 to 0.9 mm., and were conveniently reduced to 
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Fig, 4,6 E vs E curves for x(0)=-180 m/sec 2 and x(0)=.10 m/sec. 
Curve (a) is for a vq(0)=-1 00 volts, (b) for 
vc(0)=-85 volts. Notice in this curve that although 
there is a tendency to override, because of a high speed 
of entry, the system is kept within the bounded region 
by the adequate use of vc(0)=0. Curve (c) shows an 
overshoot caused by the use of vq(0)=0. 

x scale = 5xl(T 5 m per inch y scale = ,02 m/sec per i 
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0.05 or 0.04. The integral squared error was reduced by a factor 
which ranges between 20 and 300, which means that the performance of 

the whole system is noticeably improved. 

Finally, it is important to point out from the experimental data, 
that when reducing the integral squared error, not only does the 
trajectory stay inside the bounded region, but also the settling time 

is minimum for this case. 
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V . CONCLUSION 



The special case presented here has served to show the significance 
of the initial condition of a filter when connected to a system to 

e 

compensate it. In practice, other types of filters can be used, and 
a different type of feed encountered, but the analysis will be similar, 
and the effects of the initial condition should be carefully studied 
in order to obtain an optimized control. 

The bounded characteristics of this system are evidently no 
restrictions for application of the criterion here presented. It is 
generally desired to achieve minimum time of response, thus making 
the assumptions true for all cases. More sophisticated types of control 
can be used based on the results obtained, such as adjusting the 
initial condition according with the values of other parameters of the 
system if these are well known, or other convenient control that, without 
affecting cost or other factors, will provide a significant improvement 
on system performance. 

Thus, the main idea introduced by this work is that when using a 
compensator in a discontinuous manner, initial conditions of the system 
play an important role in its response characteristics, and should 
therefore be considered carefully to obtain the best possible solution. 
The specific control law used must be determined for each particular 

case . 
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APPENDIX A 

Explanation of Simulation 

S/360 CSMP is a program for the simulation of continuous systems, 

c 

hence limiting the possibility of simulating the complete discontinuous 
problem in one single program when a filter with initial condition is 
added to the diagram after time =0. This leads to the use of two 
separate programs in which, the initial conditions of the second are 
the final values of the first. 

In order to account for the initial conditions as specified by the 
problem at hand, some of the operations had to be made longer, as is 
the case of the three operations involved in solving the filter 
equations instead of just using the LEDLAG statement provided by CSMP, 
and the two separate integrations for the motor equations in place of 

k 

the CMPXPL statement. 

The program provides an application-oriented input language that 
accepts problems expressed in the form of either an analog block 
diagram or a system of ordinary differential equations. The N0S0RT 
statement in the discontinuous mode program allows the use of logical 
statements as the ones known in Fortran. 



LABEL MULTIMODE AUTOMATIC CONTROL SYSTEM 



DISCONTINUOUS OPERATION PHASE 

PARAM KM = 3432.3, ZETA = .9985, OMEGA = 
PARAM K = 55000, OMSQ = 72740.95, TAU = . 
PARAM EINIC = ( .02544, .00254) 

£ 

OYNAMIC 

NOSORT 

INPUT = 24 

IF ( E.LE.EINIC/7 ) INPUT = -24 

SOLUTION TO INSTRUMENT SERVO EQUATIONS 

C3D0T = AUX1 + AUX2 + AUX3 
C2D0T = I NTGRL ( 0 . , C3D0T ) 

COOT = I NTGRL ( 0 . ,C2D0T) 

C = I NTGRL ( 0 . , CDOT ) 

AUX1 = -2 * ZtTA * OMEGA * C2D0T 
AUX2 = -OMSQ * CDOT 
AUX3 = KM * INPUT 
EDOT = -CDOT 
E = EINIC - C 

FINISH STATEMENT TERMINATES THE INTEGRATION FOR 
SPECIFIED VALUE OF CDOT 

FINISH CDOT = -.001 

PRINTOUT OF THE INDICATED VARIABLES 

PRINT C3DQT,C2D0T,CD0T,C,E 

DATA FOR DRAW SUBROUTINE 

PREPARE E , EDOT 

TIMER FINTIM = .1, CELT = .001, PRDEL = . 



ENDJOB 



END 

STOP 



269.7 

004 



001 
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LABEL MULTIMODE AUTOMATIC CONTROL SYSTEM 



CLOSED LOOP, LINEAR OPERATION PHASE 

PARAM KM = 3432.3, ZETA = .9985, OMEGA = 269 
PARAM K = 55000, OMSQ = 72740.95, TAU = .004 
INCCN CO = -6.35E-05 

INCON C2D0T0 = -160, CDOTO = .08, VC0=( -100,- 
DYNAMIC 

EPRIME = E * K 

SOLUTION TO FILTER EQUATIONS 

VC IS THE VOLTAGE ACROSS THE CAPACITOR 
EC IS THE OUTPUT OF THE COMPENSATOR NETWORK 

VCDOT = -( 1 0/ TAU ) * VC + ( 9/TAU ) * EPRIME 
VC = INTGRLIVCO, VCDOT) 

EC = EPRIME - VC 

SOLUTION TO INSTRUMENT SERVO EQUATIONS 

C3D0T = AUX1 + AUX2 + AUX3 
C2D0T = INTGRL ( C2DOTO , C3D0T ) 

CDOT = I NTGRL (CD0T0,C2D0T) 

C = INTGRL (CO, CDOT ) 

AUX 1 = -2 * ZETA * OMEGA * C2D0T 
AUX2 = -OMSQ * CDOT 
AUX3 = KM * EC 

E2D0T =- C2D0T 
EDOT = - CDOT 
E = - C 

EVALUATION OF THE INTEGRAL SQUARED ERROR 

SQER = E * E * l.OE 10 
INSQER = INTGRH 0. , SQER) 

PRI NTOUT OF THE INDICATED VARIABLES 

PRINT E2D0T , E DOT , E , I NSQER 

DATA FOR DRAW SUBROUTINE TO OBTAIN GRAPHICAL OUTPUT 
OF E VS. EDOT 

PREPARE E , EDOT 

TIMER FINTIM = .1, DELT = .001, PRDEL = .001 

END 

STOP 

ENDJOB 



7 

85,0 ) 
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CSMP 



TIME 

0.0 

1 .00 00 E- 03 
2. 0000E-03 
3. OOOOE-03 

4 . 0000 E— 03 
5. 0000E-03 

6.0000 E— 03 
7. OOOOE-03 
8. 0000E-03 

9 . 0000 E— 03 
1 . 0000E-02 
I. 1000E-02 
1 .2000E-02 
1 .3000E-02 

1 .4000 E— 02 

1 . 5000 E— 02 
1.60C0E— 02 
1.7000E-02 
I . 8000E— 02 
1.9000E-02 

2.0000 E— 02 
2. 1000E-02 
2. 2000 E-02 
2. 3000E-02 

2.4000 E— 02 

2.5000 E— 02 
2 . 6000 E-02 
2 . 7000E— 02 

2 . 8000 E— 02 

2 . 9000 E— 02 

3.00 00 E- 02 
3. 1000 E-02 
3. 2000 E-02 
3 . 3000 E- 02 

3 . 4000 E-02 
3. 5000 E-02 
3.6000E-02 
3 . 7000E-02 
3. 8000E— 02 

3.9000 E-02 
4. 0000E-02 
4. 1000 E-02 
4.2000E-02 
4. 3000E-02 

4.4000 E— 02 
4. 5000E— 02 
4 . 6000E— 02 
4.7000E-02 

4. 8000 E-02 

4 . 9000 E— 02 
5 . 0000E-02 
5. 1000E-02 
5 .2000 E-02 
5 . 3000E-02 
5. 4000 E-02 

5 . 5000 E— 02 



E2D0T 

1.6000E 02 
9. 1178E 01 
4.4269E 01 
1.0041E 01 
-1.3587E 01 
-2.8262E 01 
-3.5644E 01 
-3.7365E 01 

— 3 . 4956 E 01 
-2.9794E 01 
-2.3063E 01 
-1.5734E 01 
-8.5570E 00 
-2.0755E 00 

3.3608E 00 
7.5678E 00 
1.0498E 01 
1.2210E 01 
1.2836E 01 
1.2560E 01 
1.1584E 01 
1.0119E 01 
8.3597E 00 
6.4805E 00 
4.6254E 00 
2. 90 58 E 00 
1.3999E 00 
1 . 5505E-0 1 
-8.0851E-01 
-1.4932E 00 
-1.9187E 00 
-2.1172E 00 
-2.1279E 00 
-1.9935E 00 
-1.7564E 00 
-1.4559E 00 

— 1 . 12 66 E 00 
-7. 9684E-01 
-4. 8874E-01 
-2. 1780E-01 

6.4451E-03 
1 . 7974E-0 1 
3.0215 E— 0 1 
3. 7712 E— 0 1 
4. 1042 E— 0 1 
4. 0928E-01 
3. 8157E-01 
3. 3511E-01 
2.772CE-01 
2. 1424E-01 
1.51 50 E— 0 1 
9. 3059E-02 
4. 1778 E— 0 2 
-5.9734E-04 
—3 . 3303 E-02 
-5.6383E-02 



EDOT 

-8. 0000 E-02 
4. 2633E-02 
1.0916E-01 
1.3536E-01 
1 .3277E-01 
1. 1117 E— 01 
7 . 8680E— 02 
4. 1770E-02 
5.3252 E-03 
-2.7227E— 02 
-5.3742E-02 
—7. 3156E— 02 
— 8 . 5263E-02 
-9. 0504 E-02 

— 8.97 64E— 02 
-8.4194E-02 

— 7. 5056 E- 02 
-6. 3605 E-02 
-5 . 0999E-02 
—3. 8234E-02 
-2.6112E-02 
-1. 5228E-02 

— 5 . 9720E- 03 
1.4516E-03 
6 . 9974E— 03 
1.07 48E- 02 
1.2881 E— 02 
1.3635E-02 
1 . 3285 E-02 
1 . 2 1 1 1 E— 02 
1 . 0385E— 02 
8 . 3496 E— 03 
6 . 2132 E— 03 
4. 1422E-03 
2.2603E-03 
6. 5040E-04 

-6.4198E-04 

— 1 . 60 27 E-03 
—2. 2429E— 03 
-2 .5926E-03 
- 2 . 6942 E-03 

— 2 . 596 8E- 03 
-2 . 35 17 E-03 

— 2 . 00 83 E— 03 
—1.61 14E- 03 
-1. 1990 E-03 
-8.0 167E-04 
-4. 4209E— 04 

— 1 . 3526 E— 04 
1. 1064E-04 
2.9331 E— 04 
4. 1509E-04 
4 . 81 83 E— 04 
5.0163 E— 04 
4. 83 86E-04 
4.3824E-04 



E 

6.35 00E-05 
5. 0477E— 05 
1.3027E-04 
2. 5 5 38 E— 04 
3.9142E-04 
5 . 146 1 E— 04 
6.1015E-04 
6.7052E-04 
6.93 86 E— 04 
6. 8248E— 04 
6.4143E— 04 
5.7737 E— 04 
4. 9756E-04 
4.09 14E-04 
3. 1855E-04 
2.3122E-04 
1.5135 E— 04 
8. 1881E-05 
2 . 45 27 E— 05 
-2.0066 E-0 5 
— 5.2157E-05 

— 7.2705 E— 05 
—8. 3158E-05 

— 8. 5262 E— 05 
—8 . 08 83 E-05 
—7. 1867E-05 

— 5 .9927 E— 05 
-4.6565E-05 
-3. 3025E— 05 
—2. 0270E-05 

— 8.98 64 E— 06 
3. 9736 E-07 
7 .6796 E— 06 
1. 2846E-05 
1 . 60 27 E— 05 
1.7458 E— 05 
1.7434E-05 
1 . 62 85 E— 05 
1. 4336E— 05 
1.1896 E-0 5 
9 . 23 36 E— 06 
6. 5737E-06 
4. 0893E— 06 
1 . 9030 E— 06 
9. 04 10 E— 08 

— 1.3147 E— 06 
-2.3127E-06 
-2.9307 E— 06 
-3.2145E-06 
-3. 22 16 E — 06 
—3. 0144E— 06 
-2.6553E-06 
-2. 2026E-C6 
-1.7073E-06 
— 1.2118 E-06 
-7. 48 85 E-07 



INSQER 

0.0 

2. 2244E-02 
9. 97 13 E-02 
4. 7640 E— 01 
1.5392E 00 
3.6205E 00 
6.8214E 00 
1.0964E 01 
1.5660E 01 
2.0433E 01 
2 .4846 E 01 
2.8583E 01 
3.1488E 01 
3.3554E 01 
3.43 84E 01 
3.5643E 01 
3.6011E 01 
3.6149E 01 
3.6179E 01 
3.6181E 01 
3.6195E 01 
3.6236E 01 
3.6298E 01 
3.6370E 01 
3.6440E 01 
3.6499E 01 
3.6542E 01 
3.6571E 01 
3.6587E 01 
3.6594E 01 
3.6596E 01 
3.6596E 01 
3.6597E 01 
3.6598E 01 
3.6600E 01 
3.6603E 01 
3.6606E 01 
3.6609E 01 
3.6611E 01 
3.6613E 01 
3.6614E 01 
3.6614E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 
3.6615E 01 



RKS 



INTEGRATION 



VCO 



0.0 



CSMP 



TIME 

0.0 

1.0000E-03 
2 • 0000E-03 

3.0000 E— 03 
4. OOOOE-03 

5 . 00 00 E— 03 
6. 0000E-03 
7. 0000E-03 

8.00 00 E-03 

9 . 0000 E— 03 

1.0000 E— 02 
1 . 1000E-02 
1.2000E-02 
1 . 3000E-02 
1.4000E-02 
1 . 5000E-02 

1.6000 E— 02 
1 .7000 E-02 
1. 8000E-02 

1.9000 E— 02 

2.0000 E-02 
2. 1000 E— 02 

2.2000 E-02 
2 . 3 000 E— 02 
2.4000E-02 
2 . 5000 E-02 
2. 6000E-02 
2 . 70 00 E-02 

2 . 8000 E-02 
2 . 9000E— 02 
3.0000 E-02 
3. 1000E-02 

3 . 2000 E— 02 
3.3000 E-02 
3.4000E-02 
3. 5000E-02 

3.6000 E— 02 
3 . 7000E-02 

3 . 8000 E— 02 
3.9000E-02 

4.0000 E— 02 
4. 1000E-02 
4. 2000E-02 
4. 3000 E— 02 
4.4000E— 02 
4. 5000E-02 
4. 6000E— 02 
4.7 000 E— 02 

4 . 8000 E— 02 

4 . 9000 E— 02 
5 . 0000E-02 
5. 1000 E-02 
5 . 2000E-02 
5 . 3000E-02 
5. 4000 E-02 
5. 5000E-02 



E2D0T 

1.6000E 02 
6.1846E 00 
-9 . 0044E— 0 1 
2.1396E 00 
4.6667E 00 
6.0623E 00 
6.5067E 00 
6. 23 51 E 00 
5.4687E 00 
4.4023E 00 
3. 1991 E 00 
1.9890E 00 
8.6843 E— 0 1 
-9.7302E-02 
-8.7046E-01 
-1.4366E 00 

— 1 . 7998 E 00 
-1.9779E 00 
— 1.9975E 00 
-1.8908E 00 
-1.6913E 00 
-1.4316E 00 
-1.1411E 00 
-8. 4496E— 01 
—5 . 63 19E-0 1 
-3. 1045E-01 
-9. 6309E-02 

7 . 426 IE-0 2 
2.002 IE— 0 1 
2. 8360E-01 
3.2878 E— 0 1 
3.4158 E— 0 1 
3.2861 E— 0 1 
2 . 9666E— 0 1 
2. 5220E-01 
2 . 01 05E-0 1 
1.48 14E— 0 1 
9.7371E-02 
5. 1630E-02 
1 . 2804E— 02 
-1. 8106E-02 
-4. 0851 E-02 
-5. 5777 E-02 
-6. 3670E-02 
-6. 5600E-02 
-6. 2789E-02 
—5 . 6496E— 0 2 
—4. 7922E-02 
—3. 8146E— 02 
-2 . 8084E-02 

— 1 . 846 1 E-0 2 
-9. 8099E-03 
-2. 47 92 E— 03 

3. 3490E-03 
7.63 26E-03 
1 . 0440E— 02 



E00T 

— 8.00 OOE— 02 

— 2 . 7655 E— 02 
-2.7320 E— 02 
—2. 6765E— 02 
-2 . 3266E-02 

— 1.781 3 E— 02 
-1. 1459E-02 
-5.0381 E-03 

8.4653E— 04 
5 . 7997E— 03 
9 . 6059 E— 03 
1 .2196E-02 
1.3614E-02 
1.39 85 E — 02 
1 . 3484E— 02 
1 . 23 13 E— 02 
1 . 0679E— 02 
8 .77 53E-03 
6 .7757 E— 03 
4 . 8224E— 03 
3.0250E-03 
1.45 99 E — 03 
1 . 7209E-04 
—8. 2050E— 04 

— 1 . 5227E-03 
-1.9566E-03 

— 2.1565 E— 03 
-2. 1638E-03 
-2.0229E-C3 

— 1 . 7776E— 03 
-1.4684E-03 

— 1.13 08E— 03 

— 7.938 9E— 04 

— 4.799 5 E-04 

— 2 . 0473 E— 04 
2. 2232E-05 
1 . 9680E— 04 
3.1924E— 04 
3.9323E-04 
4. 248 1 E— 04 
4.2148E-04 
3.9133E-04 
3.4239E-04 
2. 8212E-04 
2.17 04E— 04 
1 . 5250E— 04 
9 . 2623E— 05 
4.0271 E-05 

—2. 8222E-06 
-3.5928E-05 

— 5 . 91 38 E— 05 
-7 . 31 75E-05 
-7. 9200 E-05 
—7. 8636E— 05 
-7.3018E-05 
-6. 3864 E-05 



E 

6. 3500E— 05 
2. 1242E-05 
— 5 . 7541 E-06 
-3. 3057E-0 5 
-5.82 83 E— 0 5 
-7. 8939E-05 
-9. 3612E-05 
-1.01 84E-04 
-1. 0387E-04 
-1. 0046E-04 
-9. 2654E-05 
-8. 1652E-05 
-6 . 86 53 E- 05 

— 5 . 4773 E— 0 5 
-4. 0975E-05 
-2. 8029E-05 
-1.65O3E-05 
-6 .76 10E-06 

1.0161 E-06 
6. 8 062 E— 06 
1.0713 E-0 5 
1. 2934E-05 
1. 3726E-05 
1. 3377E-05 
1.2182E— 05 
1. 0421E-05 
8.3466E-06 
6. 1723E-06 
4.0684E-06 
2. 1612E-06 
5. 3445E-07 
—7. 6 62 5 E— 07 
-1. 7275E-06 

— 2 . 36 1 8 E— 06 
-2.7004E-06 
-2.7874E-C6 
-2. 6735E-06 
-2.4112 E-06 
-2.05 12E-06 

— 1 . 63 89 E-06 
-1. 2132E-06 
-8.0489E-07 
-4. 3679E-07 
-1. 2388E-07 

1.2587E-07 
3.1 040 E— 07 
4. 3244E-G7 
4.98 17 E-07 
5. 1608E-07 
4. 9587E— 07 

4. 47 53 E— 07 
3. 8065E-07 
3.03 86E— 07 
2 . 2445 E-07 

1 . 48 27 E— 07 
7. 9593E-08 



INSQER 

0.0 

1 .6137E-02 
1 . 7391E— 02 
2. 1798E— 02 
4. 3457E-02 
9.15 15E-02 
1 . 67 05E-01 
2 . 6365 E- 0 1 
3.7 046 E— 0 1 
4 . 7569 E— 0 1 
5.6958E-01 
6. 4601E— 01 
7 . 028 IE-01 
7.41 09E— 01 
7. 6413E-01 
7. 7611 E— 01 
7.811 2E— 0 1 
7.825 1 E— 01 
7. 8264E-01 
7. 82 83 E— 0 1 
7 . 8364E— 0 1 
7.85 07 E— 0 1 
7. 8687E-01 
7. 8873E-01 
7.9038E-01 
7 . 9167E-01 
7.9256E-01 
7. 9309E-01 
7 . 933 5 E— 0 1 
7 . 9345 E— 01 
7. 9347E-01 
7.9347E-01 
7 . 9349 E— 0 1 
7. 9353E— 01 
7 . 9360E— 0 1 
7.9368E-01 
7.9375E-01 
7. 9382 E- 01 
7.9387E-01 

7 . 93 90 E-0 1 
7 . 9 392 E-0 1 
7 . 9393E— 01 
7 . 9393E— 01 
7. 9394E-01 
7. 9394E-01 
7.9394E-01 

7.93 94E-0 1 
7.9394E-01 
7.9394E-01 

7.93 94E-0 1 

7 . 93 95 E-0 1 
7. 93 95 E— 0 1 
7. 9395E— 01 

7.93 95 E— 0 1 

7 . 93 95E-01 

7 . 93 95 E-0 1 
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CSMP 



TIME 

0.0 

1.00 00 E-03 

2.00 00 E— 03 

3 . 0000 E— 03 

4.0000 E— 03 
5.0000E-03 
6. 0000E-03 
7. 0000 E-03 

8 . 00 00 E— 03 
9. 0000 E-03 

1.00 00 E— 02 

1 . 1 000 E— 02 
1. 2000E-02 

1.30 00 E— 02 

1 . 4000 E— 02 

1 . 5000 E— 02 
1 .60 00 E— 02 
1.7000 E-02 

1 . 8000 E— 02 

1 . 9000 E— 02 
2. OOOOE-02 
2. 1000 E-02 
2 . 2000E-02 
2 . 3000 E-02 
2.4000E-02 
2. 5000 E-02 
2. 6000 E-02 
2.7 000E— 02 
2. 8000E-02 
2.90 00 E— 02 

3 . 0000 E— 02 
3 . 1000E-02 
3 . 2000 E-02 

3.30 00 E— 02 
3.4000E-02 

3 . 5000 E-02 
3. 6000E— 02 
3. 7000 E-02 

3 . 8000 E— 02 

3 . 9000 E— 02 
4. 0000E-02 
4. 1000E-02 
4. 2000 E— 02 
4.3 000 E— 02 

4 . 4000 E— 02 

4.5000 E— 02 
4 . 6000 E— 02 
4. 7000 E-02 

4.8000 E— 02 

4. 9000 E— 02 
5 . 0000E-02 

5.1000 E-02 
5. 2000 E-02 
5. 3000 E-02 
5. 4000 E— 02 

5 . 5000 E— 02 



E2D0T 

1.6000E 02 
1.8934E 01 
5.8752E 00 
3.3248E 00 
1.9287E 00 
9. 1371E-01 
1. 8405E-01 

— 3. 049 IE— 01 
—5 . 95 00 E— 0 1 
-7. 2719E-01 
—7. 4023E— 01 
-6. 6937E-01 
-5.4539E-01 
-3.9403E— 01 
—2. 3577E— 01 
-8. 5958E-02 

4 . 4837 E— 0 2 
1. 5027E-01 
2. 2757E-01 
2.7678 E— 01 
3.0006E-01 
3.0101E-01 

2. 8404E-01 
2. 5386E-01 
2. 15 1 IE— 0 1 
1.7199E-01 
1. 2812E-C1 
8.6381E-02 
4. 89 03 E— 0 2 
1.7083E-02 
-8. 3505E-03 
-2. 7237E-02 
-3. 9869E-02 
—4. 6872E— 02 
-4.9094E-02 
-4 .7496 E— 02 
-4. 3070E-02 
-3. 6762E-02 
-2. 9426E-02 
-2. 1787E-02 

— 1 . 44 24E— 0 2 
-7. 7627 E-0 3 
-2. 0877E-03 

2 . 4484 E— 0 3 
5. 8030E-03 
8.02 13 E— 0 3 
9. 2138E-03 
9. 5335E-C3 
9. 1565 E— 03 
8. 2650E-03 
7. 0338 E— 03 
5. 62 07 E— 03 
4. 1595E-03 
2. 7571 E— 03 
1.4923E-03 
4. 1684E-04 



EDOT 
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